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A new form of a binary Darboux transformation is used to generate analytical solutions of a 
nonlinear Liouville-von Neumann equation. General theory is illustrated by explicit examples. 
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I. INTRODUCTION 

56 ! 

Q\ • Nonlinear operator equations one encounters in quantum optics and quantum field theory are typically solved 
] by techniques which are either perturbative or semiclassical (cf. QB]). The situation is caused by the fact that 
. analytic methods of dealing with "non-Abelian" nonlinearites are still at a rather preliminary stage of development. 
An important step towards more efficient analytical techniques is associated with the notion of an inverse spectral 
transformation. The use of the method in the contex of matrix equations can be found in ||-||] where an analytical 
treatment of Maxwell-Bloch equations is given. In application to the Maxwell-Bloch system describing three-level 
CNj | atoms interacting with light || one makes use of a degenerate Zakharov-Shabat spectral problem with reduction 
constraints |Q . The same problem is used in the context of the complex modified Korteweg-de Vries (MKdV) equation 
for a slowly varying envelope of electromagnetic field in an optical fiber Q . 

A technical complication occurs if a solution obtained by an inverse method should additionally satisfy some 
' constraint. For example, it is often essential to guarantee that the solution one gets is Hermitian or positive. Difficultes 
, of this kind were one of the motivations for the development of new Darboux-type operator techniques of solving non- 
Abelian equations. Particularly useful turned out to be the method of elementary and binary Darboux transformations 
introduced by one of us |p|-pT| . These particular versions of the Darboux transformations are more primitive than 
the ordinary ones in the sense that the latter can be obtained by their composition. The binary transformation, 
a result of an application of two mutually conjugated elementary Darboux transformations one after another, was 
successfully applied to a three-state Maxwell-Bloch system with degeneracy in |^] , and various multisoliton solutions, 
4 ' including the well known 27r-pulse and breathers, were found. 

Qh, In this article we apply a generalization of this technique to a new type of nonlinear operator equation. The 
i— > nonlinear Liouville-von Neumann equation we will discuss is the simplest nontrivial example of a Lie-Nambu dynamics 
of a density matrix and occurs naturally in certain version of nonlinear quantum mechanics. To begin with, let us 
recall that the well known Liouville-von Neumann equation (LvNE) 



ip=[H,p], (1) 

where H is a Hamiltonian operator, p a density matrix, and the dot denotes the time derivative, is linear. In 
Hartree-type theories one considers more general, nonlinear equations of the form 

ip=[H(p),p], (2) 

where H(p) is a nonlinear Hamiltonian operator. For time- independent Hamiltonians H yp(t)) = H (p(0)) the formal 
solutions are exponential 

p(t)=exp[-iH(p(Q))t]p(0)exp[iH(p(0))t\, (3) 

Both kinds of nonlinear LvNE's can be written in either Lie-Poisson JlJ,|lJ] or Lie-Nambu forms fl5|-|l9||. The Lie- 
Nambu version involves a 3-bracket and the LvNE's can be written as 

ip a = {p a ,H f ,S} = {p a) H f } (4) 

where {•, •} := {•, •, S} is a Lie-Poisson bracket. Here p a :— pAA'( a , a ') are components of p in some basis, A and A' 
are discrete (say, spinor) indices and a, a' the continuous ones. Hf is a Hamiltonian function and S a functional that 
can be identified with the 2-entropy of Daroczy and Tsallis p2| , i.e. 2S = Tr(p 2 ) =|| p || 2 is the Hilbert-Schmidt 
squared norm of p. 

An extension from a Lie-Poisson bracket to a 3-bracket led Nambu to a generalization of classical Hamiltonian 
dynamics |23|] . The 3-bracket equation (||) naturally leads to the question of possible Nambu-type extensions of the 
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Lie-Poisson dynamics of density matrices. An interesting class of such generalizations occurs if one keeps Hf(p) linear 
in p but S is a function of other Daroczy-Tsallis entropies. Such Nambu-type equations are rather unusual from the 
point of view of generalized Nambu-Poisson theories (24| . The peculiarity is that although the 3-bracket itself does 
not satisfy the so-called fundamental identity, typically regarded as a Nambu analogue of the Jacobi identity, the 2- 
bracket defined via {•,•}#/ := {-,Hf, •} does satisfy the ordinary Jacobi identity if Hf(p) is a linear functional of p Jl8| . 
It follows that restricting Hf to linear functionals one effectively uses the Nambu-type structure as an intermediate 
step which allows one to introduce a new Poisson structure, and now 5* plays a role of a Hamiltonian function. The 
standard way of introducing interactions, i.e. by modifying Hf, looks from this perspective as a modification of the 
Poisson structure while keeping the Hamiltonian S fixed. A simultaneous change of the two Nambu-type generators, 
Hf and S, can be regarded as a change of the Hamiltonian function S, accompanied by a modification of the Poisson 
structure defined in terms of Hf. In this respect the Nambu-type bi-Hamiltonian dynamics with linear Hf has a logical 
structure analogous to this of general relativity. The choice of linear Hf and generalized 5* can be also motivated by 
difficulties with probability interpretation of generalized observables since there is no physically natural definition of 
spectrum of nonlinear operators |ll|[35|,[5(| . The 3-bracket structure can be shown to be a a particular case of a still 
more general (2/c + l)-bracket one that, for 2k + 1 > 3, always vanishes on pure states and therefore is invisible at the 
level of the Schrodinger dynamics p8[ . 

The nonlinear LvNE corresponding to S — Tr(p™)/n, 

ip=[H,p n -% (5) 

was introduced in Jl6| . General properties of such equations were discussed in | jl7| and ]l8| . It was shown, in particular, 
that spectra of their Hcrmitian Hilbcrt-Schmidt solutions are time-independent. This opens a possibility of a density 
matrix interpretation of the solutions. Let us note that for p 2 = p (pure states) the equations reduce to the linear 
LvNE and, therefore, the pure state dynamics is indistinguishable from the ordinary linear Schrodinger one. One of 
the problems that still remained open was how to solve such nonlinear equations. There exist formal solutions given 
in a form of a series, but the question of convergence of such a series was not investigated. 

The aim of this paper is to describe an algebraic method that leads to solutions of a nonlinear LvNE which reduces 
to (||) with n = 3 (3-entropy equation) in special cases. The equation we shall study is 

ip = [H, p 2 } + ip'H + iHp' (6) 

where the prime denotes a derivative with respect to some additional parameter r. We will generate the solutions 
from a Lax pair with the help of a binary Darboux transformation. To avoid technicalities we will generally assume 
that the Hamiltonian H and other operators are finite-dimensional matrices, but the transformation works in a much 
more general setting (see the example of the harmonic oscillator) and its application to general infinite-dimensional 
systems is a subject of current study. 



II. LAX PAIR AND ITS DARBOUX COVARIANCE 



The technique of Darboux-type transformations is perhaps the most powerful analytical method of solving differ- 
ential equations. Although it was developed mainly in the context of nonlinear equations, it is implicitly used also in 
standard textbook quantum mechanics under the name of the creation-annihilation operator method. The method of 
creation operators is simultaneously a good illustration of the way the Darboux technique works. In short, to use the 
method one has to begin with an initial solution which is found by other means (a "ground state"). Then one has 
to find a "creation operator" and the Darboux transformation is a systematic procedure that allows one to do it. In 
linear cases once we have these two elements, we are able to generate an entire Hilbert space of solutions. In nonlinear 
cases the spaces of solutions are bigger and therefore a given "ground state" and a "creation operator" may generate 
only a subset of this space. It is mainly for this reason that much effort was devoted to finding different generalizations 
of Darboux transformations (cf. [fl2"[ ) . The method we will use was devised for non-commutative equations such as 
Heisenberg equations of motion. The construction given in |lC| , |ll] ] led to a transformation more general than the one 
we use and its derivation from elementary transformations is somewhat tedious. However, once one has our explicit 
form, one can check by a straightforward calculation that the binary transformation indeed maps one solution into 
another. To make this paper self-contained we give the explicit proof in the Appendix. 

Consider the following pair of Zakharov-Shabat equations 

icp(p,)' = (U- (iH)cp(p,) =: Z&Qi) (7) 
i(p(p,) = (UH + HU - pH 2 )tp{p) (8) 
= \{IJ 2 -Zl)^p) (9) 
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where U and H are Hermitian matrices the dot and prime denote, respectively, derivatives with respect to time t and 
some auxiliary parameter r, and p, is complex. The solution ip(p) is also in general a matrix. We assume that H is t 
and r-independent and U — U(t,r). The compatibility condition for (^), (JsJ) is 

ill = [H,U 2 ]+iU'H + iHU', (10) 

and therefore the above pair is the Lax pair for (Q). We will stick to the notation with U instead of p since non- 
Hermitian and non-positive solutions are also of some interest and p will be reserved for density matrices. 
We will need two additional conjugated problems 

- iip(X)' = i/j(X)(U - XH) (11) 
-iip(X) =ip(X)(UH + HU - XH 2 ) (12) 



-%x{v)'=x{v)(JI-vH) 
-ijdy) = x{v){UH + HU- vH 2 ) 

each of them playing a role of a Lax pair for (^) . 

Consider for the moment the following general Zakharov-Shabat problems 

idip(p) = (V - pJ)ip(p) 
-idip(X) = ip(X)(V - XJ) 
-idx{v) = x{v){V - vJ) 



(13) 
(14) 



(15) 
(16) 
(17) 



where d denotes a derivative with respect to some parameter. We will take the binary transformation in the form 

</>[l](A,/i,i/) 



= V(A) 



v — fi 
A — fi 

v — fj, 

X — ji 



P 



(18) 
(19) 



where p is a constant projector (dp — 0) and the inverse means an inverse in the p- invariant subspace: [pxp)~ 1 pxp = 
pxp(pxp)^ 1 = p. The operator P defined by ( fl9|) is idempotent (P 2 — P) but in general non-Hermitian. P satisfies 
the nonlinear master equation j37[ 



idP = (V - pJ)P - P(V - vJ) + (/i - v)PJP 
The binary transformation implies the following transformation of the potential 

V[l](n,v)=V + (n-v)[P,J}. 
Applying this general result to V = U, J = H we get 

U[l]{fru) = U + (ji-v)\P,H\. 
The second triple of equations we have started with corresponds to V = UH + HU and J = H 2 . In this case 

V[l](ji, v) = U[1](m, v)H + HU[l]Qi, v). 



(20) 
(21) 
(22) 
(23) 



This means that ( p2[ ) guarantees simultaneous covariance of the Lax pairs under the binary transformation (18). 

Another important feature of the binary transformation is the fact that for v = p and px(P-) — the 
Hermiticity of the potential is Darboux-covariant, i.e 

U[l](p,tf = U[l](ji,ii) (24) 

if [/t = U. 

Using ( pp[ ) one can show by a straightforward calculation (see the Appendix) that the binary transformed 

^[l]=V>(l-^V) (25) 
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indeed satisfies 

-im' =m{u[i]-\H) (26) 

-#[1] = iP[l]{U[l]H + HU[1] - \H 2 ) (27) 

with U[l] = U[l](p,ft) and, therefore, 

iU[l] = [H, U[l] 2 ] + iU[l]'H + iHU[l}'. (28) 

Subsequent iterations of the Darboux transformation generate further solutions. Starting with a Hermitian solution 
U we obtain an infinite sequence of Hermitian solutions U[l], U[2],... satisfying Tr U —Tt U[l] =Tr t/[2] ... 



III. COVARIANCE OF THE CONSTRAINT U'H + HU' = 

In order to generate solutions of ill = [H, U 2 ] one has to maintain the constraint U'H + HU' = U[1]'H + HU[1]' = 
U[2]'H + HU[2]' = , . . = 0. Starting with stationary solutions 

i<p(n)' = zip(n) (29) 

one finds that U' = implies U[l]' = 0. An alternative approach can be applied to Hamiltonians of the Dirac type 

H = p ■ a + m0, (30) 

which satisfy H 2 — E(p) 2 l and therefore imply 

U[1]'H + HU[1}' = U'H + HU' + (m - P)[P', H 2 ] 

= U'H + HU' (31) 

which makes the constraint hereditary. 

In general, using (|2C|), one finds that the constraint is hereditary if 

[P X (U - fiH)P - P(U -vH)P ±1 H 2 }=0, (32) 

where Pi = 1 P. 



IV. PARTICULAR CASES 



In this section we shall discuss properties of solutions corresponding to several choices of the initial U. 
For some applications one can restrict the general form (|l8|) by choosing 



0. 



v — fi and x(/x) = y(/Lt)t. Denoting the first column in by \<p) one finds that 

U[i\0hf£) = U+0i-fi)[P ) H\ 

where 



\ 


/ <pl 


0.. 


• \ 






0.. 


. 


, ^(m) = 






/ 




0.. 


■ °/ 



(33) 



P = 



{<P\<P) 



(34) 



(35) 



The transition from the column solution \ip) to ( ^3| ) is a useful trick that allows one to consider expressions such as 
pipp which otherwise would not make any sense. 
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A. U 2 = U, V = 



This is an interesting case since U = U(t) is a solution of the ordinary linear LvNE: 

U(t) = exp [ - iHt] P(0) exp [iHt] . (36) 
Let us take a solution stationary with respect to r: 

= (U- tiH)\ip{n)) = (37) 



and define 



The Lax pair is now 



with the solution 



\<p) =exp[iHt]\cp). (38) 

z^) = {U{0)-tiH)\(p) (39) 
#> = jz(*p-*j!)l£> (40) 



\<p(t,T)) = e -« m +"^ \<p(0,Q)) (41) 

where «t iT = ^z M (l — z M )i + z M r. The projector ( [35| ) is r-independent and satifies the linear LvNE. This implies 
that U[l] satisfies the same linear equation as U. The following lemmas explain the origin of this effect. Consider the 
general P defined by @ and V[l] = V[l]{ji,v). 

Lemma 1. dP = implies 



V [l] 2 = V 2 + (fi - v) [P(JV + VJ - vJ 2 )P± 

- Pj_(JV + VJ - /iJ 2 )p). (42) 

Proof: ( p0|) implies 

P(y - vj) -{V - ^J)P = (jj, - v)PJP (43) 

and 

(p - v){PJPJ + JPJP - JPJ) 

= [P,V]J+ J[P,V] +fiJ 2 P- vPJ 2 . (44) 



The latter formula leads directly to (|42]).D 

Lemma 2. Assume P' = and P is given by @ with V = PC/ + UH, J = H 2 . Then 

P[l] 2 = P 2 - (/i - i/)iP. (45) 

Proof : 

-iP ± P = -Pi (PP + HU - [lH 2 )P (46) 
-iPPjL = P(UH + HU - vH 2 )P ± (47) 

and therefore 

P[l] 2 = P 2 - «(// - ^)(PlP + PPj_) 
= P 2 -i{fi-u)P. 

□ 

An immediate consequence of Lemma 2 is 

Lemma 3. Assume P' = and P 2 = U. Then P[l] 2 = P[l] if and only if iP = [H,P] i.e. P satisfies the linear 
LvNE. 
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B. Case U 2 - aU / const-1, [U 2 - aU, H] = 0,U' = 



Let a be a real number. [U 2 — aU, H ] = implies 

U(t) = exp [ - iaHt] [7(0) exp [iaHt] . (48) 

Repeating the steps from the previous subsection we obtain the Lax pair 

z ll \(p) = {U{Q)- l iH)\(p) (49) 
#) = ±(A Q + a^-^)|$ (50) 

where A a = U(0) 2 — aU(0). The projector P is r-independent but possesses a nontrivial ^-dependence which follows 
from the fact that fi — p, ^ 0. Define the function 

F a (t) = (p(0, 0)| exp (i^-=^A t) |p(0, 0)) (51) 

satisfying 

(^(t,T)|y>(t,r)) =exp[i(a t , T -at, T )]F a (t), (52) 
where at, r = ^z^(a — + z M r. We find finally 

U[l](t) = e-^ff/^ + ^-^F^t)- 1 (53) 

xe-* A « ( [|^(0, 0)) (<p(0, 0)|, ff] e* A " 4 ) e Mfft 

=: e- laHt U int {t)e laHt . (54) 

Let us note that what makes f/[l](t) nontrivial is essentially the presence of F(t) in the denominator. It is precisely 
this property of the binary Darboux transformation that is responsible for the soliton solutions in the Maxwell-Bloch 
case ||. 

V. EXAMPLES 

We shall now demonstrate on explicit examples how the method works. We will concentrate on the first Darboux 
transformation U[l]. Further iterations, U[2], . . ., U[n], are also interesting and their relation to U[l] is similar to this 
between solitons and multi-solitons. The problem will be discussed in a forthcoming paper. 

A. 3 x 3 matrix Hamiltonian, a = 1 

Consider the Hamiltonian 

/0 1 \ 

H = 1 (55) 
V° A/ 

and take /i — i (for real [i the binary transformation is trivial). We begin with 

\+ 2 § 

C/(0) = ( o I - ^ I ■ (56) 



| 



£7(0) does not commute with H but 



1 / 1 

U{0) 2 - 17(0) = U{tf - U(t) = — I 1 | (57) 

0-1 
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does. The eigenvalues of U(0) — iH are z± = (l±i\/2)/2 and z_ has degeneracy 2. The two orthonormal eigenvectors 
corresponding to z_ are 




\<pi)= MM ■ M> = — Ml- <■>*') 

Taking 

|^(0,0)) = i(|^i) + |^2)) (59) 

we get F(t) = cosh(t/2) and the internal part defined by (|54| ) is given explicitly by 

U int (t) 

( 1+V2 _ V2 n -l-l \ 
2 l+e' U _ 2\/2cosh(t/2) 
l-\/2 x 1 



2 T+e 1 " 2cosh(t/2) 

\ 2\/2cosh(t/2) 2cosh(t/2) 2 / 



(60) 



One can check by an explicit calculation that ( p4|) with ([35|) and ( p0|) is a Hermitian solution of 

iU[l] = [H,U[1] 2 }. (61) 



Let us note that the solution ( |60| ) corresponds to an initial condition t/[l](0) which is different from U(0) and is 
no longer block-diagonal in the basis block-diagonalizing H. This is a consequence of the fact that P is not block 
diagonal, a fact that explains the importance of the degeneracy condition for (had we chosen z+ we would have 
obtained a (2 x 2) © 1 block-diagonal P). The eigenvalues of f [1](£) are nevertheless the same as those of U(0). This 
follows immediately from the i-independence of spectrum of [/[l](f) and the fact that ?7[l](f) tends asymptotically 
to U(t) for t — > +oo. As a consequence Z7[l](f) is neither normalized (Tr{/[1] ^ 1) nor positive and hence cannot be 
regarded as a density matrix. It is, however, very easy to obtain a density matrix solution once we know U[l](t). The 
problem reduces to generating a new solution whose spectrum is shifted with respect to the original one by a number. 
This can be accomplished by a gauge transformation. Indeed, 

U[l] = e ~2i\Ht fjj^ + xi) e 2lXHt (62) 

is also a solution of (fn]), and its spectrum is shifted by A with respect to this of U[l]. Such positive solutions can 
be regarded as non-normalized density matrices and are sufficient for a well defined probability interpretation of the 
theory. Let us finally note that the fact that spectrum of a Hermitian solution is conserved by the dynamics is not 
accidental but follows from general properties of Lie-Nambu equations jl7j . 

B. 3 x 3 Hamiltonian with equally-spaced spectrum 

Consider the Hamiltonian (fc, m € R) 

(k + rn —m \ 
-to k + m , (63) 
k + m J 

whose eigenvalues are k, k + m, k + 2to, and take (i — i. We begin with a non-normalized density matrix 



i(a + \/46 + a 2 ) 

P(0) = \ 1 , (64) 

c 



satisfying 
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/ 6 \ 

p(0) 2 - ap(0) = p{tf - ap(t) =06 

\0 c(c-(i) / 



(65) 



Eigenvalues of p(0) — iH are Zq = c — i(k + m), z± = \{a± y/a 2 + 4(6 — to 2 )) — i(k + to). We need this spectrum 
to satisfy a degeneracy condition: z = z + or z = 2_ with c real and non-negative. Positivity of p(0) requires also 
that a > 0, a — \jAb + a 2 > 0. We will require that 6^0 (otherwise we will not get a nontrivial p[l]) so that the 
parameters finally satisfy < Am 2 < a 2 + 46 < a 2 . Let us note that c(c — a) = 6 — to 2 independently of the choice of 
sign in the degeneracy condition z = z±. 

Denote by \k + m) the joint eigenstate of H (with eigenvalue k + m) and p(0) — iH (with eigenvalue c — i(k + to)); 
the corresponding projector is Pfe +m = |fc + m)(k + m\. Let I3 = |fc)(fc| + \k + to) (k + m\ + \k + 2m)(k + 2m|, where 
the three projectors project on eigenstates of H. We can write 

A a = bl 3 -m 2 P k+m . (66) 

The two eigenstates corresponding to the degenerate eigenvalue c — i(k + to) are othogonal. One of them is simply 
\<pi) = \k + to); the other one is ^2) = <f>k\k) + 4>k+2m\k + 2m), where the explicit form of cpj is for the moment 
irrelevant. Now take \<p{0,0)) = A\<pi) + B\ip 2 ), \A\ 2 + \B\ 2 = 1. We find 

F a (t) = e- 2bt (l + (e 2 ™ 2 * - 1)|A| 2 ) . (67) 

The above formulas can be also written as 

H= (k + n)\k + n)(k + n\ (68) 

n— 0,m,2m 

p(0) = ^(\k)(k\ + \k + 2m){k + 2m\j + c\k + m) (k + m\ 



The solution is 



where \k + j, t) = e~ iajt \k + j). 



-^\ / 4b + a 2 (\k + 2m)(k\ + \k){k + 2m|) . 



(69) 



p(t) + 22TO(l + ( e 2m2 *-l)|A| 2 ) 



(1 



x IBI 



'->k+2m 



■<t>k)(<l>. 



fe+2 m 



<j> k )\k,t)(k + 2m,t\ 



+ -j=e m2t AB(4> k+2m + <f> k )\k, t) (k + to, t\ 

+ -^=e m2t AB(0 k+2m - 4> k )\k + 2m, t)(k + m,t\ 



-H.c.) 



(70) 



C. 1-dimentional harmonic oscillator 



We begin with the Hamiltonian 

00 

H = Y,*»>{\+n)\\+n){\+n\- (71) 

n=0 

One can directly apply the construction from the above example. We have to choose some three-dimensional subspace 
which defines p(0). Put k = \ + 1 (I, to € N), and p — i/(fuv). The solution is 
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p[l](t) = p(t) + 2im{\ + (e Wt - 1)|A| 2 ) 

x (\B\ 2 ((t> k+2m + 4>k)(<l>k+2m - (t>k)\k,t)(k + 2m,t\ 

+ -Le ujm2t AB{0 k+2m - 4> k )\k + 2m, t)(k + to, t\ 
v2 



-H 



.c.) (72) 

p[l) has interesting asymptotic properties. Assume A ^ 0. For i » ~ p(t) which suggests that the nonlinear 

effect is transient. However, for t <C 

p[l](t) w + 2im((0 fc+2m + 4> k ){4> k+2m - 0fc)|fe, t) (k + 2m, t\ - H.c.) . (73) 

It follows that the asymptotic dynamics of /o[l](f) is linear but around t — some sort of "phase transition" occurs, 
and the result of this transition is stable. Let us also note that the linear evolution is determined by exp(-iaHt) with 
\a\ > 2m and m E N . The choice of a is related to the initial condition. We obtain, therefore, an effective nonlinear 
modification of frequency of the oscillator. 

Let us finally make <j>j explicit. Assume I = 0, m = 1, a = 5, b = —4, z = z + (i.e. c = (5 + V5)/2), A = B = 1/^/2. 
Now 

P(0) = §(|*N*I + IfXfl) + ^y^liXfl - |(l§Xil + IIXll) (74) 
Wi) = If) (75) 



i->=-V^W^>- (76) 

Trp(0) = (15 + V5)/2 and the eigenvalues of p(0) are 4, 1, and (5 + V5)/2. 



D. Linear equation with nonlinear perturbation 



Assume ip e = e[H,p 2 } and define 



Then 



p = cxp{—iHt]p e exp[iHt]. 



ip=[H,p\+e[H,p 2 ]. 



(77) 



(78) 



This is a Nambu-type equation obtained by taking a linear Hamiltonian function Hf = Tr (Hp) and S = Tr (p 2 )/2 + 
eTr(p 3 )/3. Average energy is, by definition (H) = Tr (Hp) /Tr p. 

Returning to the example of the harmonic oscillator we proceed as before but now we choose p = i/(ehuj). The 
solution becomes 



Pirn) 



-i(l+ae)Ht 



p(0) + 2im{\ + (e 



2eujm t 



: (\B\ 2 (4> k+2m + 4>k)(4>k+2 m - 4>k)\k){k + 2m\ 



+ 



1 
1 



l AB(l 



fe+2 m 



■<j) k )\k)(k + m\ 



-H 



t AB(4> h+ - lm -<i) k )\k + 2m)(k + m\ 



J(l+ae)Ht 



(79) 
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The asymptotic dynamics is again linear and the frequency shift is Ato = aeu>. Let us note that according to the 
definition of (H) the eigenvalues of energy should be assumed to take values hu(l/2 + ri) and not (l + ae)tiu)(l/2 + n). 
This point is essential for the probability interpretation of such a nonlinear theory. 



E. Homogeneous modification of the equation 



The equation we have solved is non- homogeneous which implies that p i— ► const • p is not a symmetry transformation. 
This fact makes it necessary to work with non-normalized density matrices. In order to obtain a homogeneous equation 
one can utilize the fact that Tr (p n ) is time-independent (as a Casimir invariant). Define C(p) = [Tr p/Ti (p 3 )] 1 ^ 2 and 
consider 

ip = C(p)[H,p 2 }. (80) 

The equation is 1-homogeneous in p and its solutions can be obtained by the substitution 1 1— > C(p)t in the correspond- 
ing formulas given above. The multiplication of p by constants is a symmetry operation so that we can easily produce 
solutions satisfying Tr p = 1. To get the equation from the Nambu-type formalism one takes S(p) = |[Tr pTr (p 3 )} 1 ^ 2 - 



F. Two spin-1/2 particles 



The above Nambu-type formalism implies that spectra of Hermitian solutions arc time- independent. In particular, 
assuming that the nonlinear dynamics is defined for a two-particle system, the corresponding two-particle density 
matrix has time-independent eigenvalues. When it comes to reduced density matrices of the one-particle subsystems 
the situation is less simple. Assume the two-particle system is described by the Hamiltonian 



H = Hi 



Ho 



(81) 



On the one hand it is clear that traces of the reduced density matrices are time-independent. On the other hand, it 
can be shown |[L8f that 



ij t Tr x (( Tr 2 pf) = 2 Tr x ([ Tr 2 (p 2 ), Tr 2 (p)]H, 



(82) 



where Tr k — 1, 2 are partial traces. For p 2 ^ p the RHS of (|82| ) dos not in general vanish and this means that the 
eigenvalues of the reduced density matrix Tr 2 p can be time-dependent. What is interesting the average energies of 
the subsystems do not change as both Tr Hi £g> lp and Tr 1 ® H 2 p are separately conserved. It follows that although 
the two subsystems do not exchange average energy, they nevertheless exhibit some kind of collective behavior. Since 
it is difficult to investigate the effect from a general perspective, it may be instructive to consider an explicit example 
of a two-particle system whose density matrix can be explicitly calculated by the Darboux technique. 
Consider two spin-1/2 particles described by the Hamiltonian 



H 



(83) 



To make the example concrete assume that |b| = 1 and \a\ — 2. We will start with the non-normalized density matrix 



/ 5 + y/7 



P(0) 



which is written in such a basis that 







5- V7 
5 + VT5 

V 5 



\ 





(84) 



15/ 



H = 2a x ® 1 + 1 ® <r z 



12 
2 10 
0-12 
2 -1 



(85) 



Take a = 5. We find 
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A 5 = p(0) 2 - 5p(0) 



9 

1 / 9 

2 | 5 
5 



(86) 



so that [A5, H] = 0. Taking \i = % we find that p(0) — iH has eigenvalues z\ 
where Z\ has degeneracy 2. The two eigenvectors corresponding to Z\ are 



(1+0/2, z 2 = (l+3*)/2, z 3 = (l-5t)/2, 



4\/2 



15 



/ -3 + iV7\ 



4\/2 



(87) 



Assuming 



we obtain 



lv(o,o)> = -L(| Vl ) + |^)). 



^ 5 (*) = ^(e 5 ' 



and p[l](t) = exp[— 5ifft]pint(t) exp[5ii?i] where 



Pint(t) = 



/ 



5 - V7tanh2t 


13i~3V7-VTE+iVW5 

8cosh2t 
7i+3\/7-3\/T5-i\/l05 

8cosh2t 





5 + \/7tanh2t 

-15t+\/7-v / 15+»\/l05 
8 cosh It 
V7+VTE 
2cosh2t 



-13i-3\/7- x/T5- 



8cosh2t 
15»+-y/7-\/i5-i\/i05 
8cosh2t 



(89) 



^L05 -7i+3\/7-3\/T5+i\/l05 
8 cosh 2t 
V^+y/lS 
2 cosh 2t 

15tanh2t 
5-v / 15tanh2t J 



(90) 



Eigenvalues p±(k), k = 1, 2, of (normalized) reduced density matrices of the fc-th subsystems are 



1 , VTE-y/7 



20 



tanh 2t 



P±(l) = 5 ± 



, . 1 ^26 + 2x/l05 

p±(2) = - ± . 

w 2 40 cosh 2t 



(91) 
(92) 



In order to check that (|82| ) is indeed satisfied one has to use non-normalized density matrices (since the equation is 

non-homogeneous) i.e. put Tr i ((Tr 2/9 ) 2 ) = 100(p+(l) 2 +p-(l) 2 ). Aver age energies of both subsystems are for 

any t, which also agrees with general theorems. 

The above collective phenomenon is typical of higher-entropy dynamics and does not occur in Hartree-type equations 
p9| , a fact that explicitly shows that the Nambu-type dynamics exhibits properties essentially different from those 
discussed in the context of completely positive nonlinear maps in [ p0[ . This interesting problem requires further 
studies and is beyond the scope of nonrelativistic theory. 



VI. CONCLUSIONS 



We have proposed an algebraic technique of solving a nonlinear operator equation. The equation we have discussed 
can be regarded as a Heisenberg-picture equation of motion for an operator U, since writing it in the form 

ill = [H, U 2 } = [HU + UH, U] (93) 

one obtains a nonlinear Heisenberg equation with the time-dependent Hamiltonian H(U) = —HU — UH. The choice 
of non-Hermitian U (typical of the binary transformation with v ^ p) leads to non-Hermitian H, a fact that may be 
of interest for a theory of open systems. 

Restricting the initial solution U to projectors (U 2 = U) we have shown that there exists a linear orbit of the 
Darboux-transformation (t/[l] 2 = U[l] and, hence, U[l] is a solution of the linear LvNE). This shows incidentally 
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that the binary transformation can be used to generate solutions of the ordinary linear LvNE, a property that may 
find applications in other contexts. 

Looking more closely at the origin of the simultaneous covariance of both equations constituting the Lax pair, 
one can immediately write other Lax pairs whose compatibility conditions provide new nonlinear Darboux-integrable 
operator equations. For example, taking the second equation with V = H 2 U + HUH + UH 2 , J = H 3 , and assuming 
the constraint U' = one obtains the compatibility condition 



iU = [H 2 U + HUH + UH 2 , U]. 



(94) 



This highly non-Abelian nonlinear Liouville-von Neumann (or Heisenberg) equation can be solved by the binary 
Darboux transformation in a way similar to this described above. 
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VII. APPENDIX: PROOF OF DARBOUX COVARIANCE 

We will show that ( |i~9| ) satisfies 



with V[l] given by @): 



idi>[l]=m(V[l]-\J) 



-idijj[l}(\,^v) 
= V(A)(V-AJ) 



v — fl 

A — /! 



P 



v — \l 

A — fi 



V>(A) 



(V - nJ)P - P(V - vJ) + (fj, - v)PJP 



v — fl 



A — /! 

- (V - \J)P + (V - fj,J)P - (A - v)PJ 
(/x - v)PJP 



v — fl 



A — /! 

(A - fj,) JP - (A - v)PJ + (fj, - v)PJP 
^[l](A,/i^)(V-[l]-AJ). 



(95) 



[1] J. H. Eberly, Am. J. Phys. 40, 1374 (1972). 

[2] S. Stenholm, Phys. Reports 6, 1 (1973). 

[3] J. C. Eilbeck, J. Phys. A: Gen. Phys. 5, 1355 (1972). 

[4] M. J. Ablowitz et al, Phys. Rev. Lett. 31, 125 (1973). 

[5] I. R. Gabitov, V. E. Zakharov, and A. V. Manakov, Sov. Phys. JETP 86, 1204 (1984). 

[6] L. A. Bolshov, V. V. Lichansky, and M. I. Persiyantsev, Sov. Phys. JETP 84, 903 (1983). 

[7] S. V. Manakov, Sov. Phys. JETP 67, 543 (1974). 

[8] N. Sasa and J. Satsuma, J. Phys. Soc. Jap. 60, 409 (1991). 



12 



[9] S. B. Leble and N. V. Ustinov, in International Symposium on Nonlinear Theory and its Applications (NOLTA '93), p. 
547 (Hawaii, 1993). 

A. A. Zaitsev and S. B. Leble, Rep. Math. Phys. 39, 177 (1997). 
S. B. Leble, Computers Math. Applic. 35, 73 (1998). 

V. B. Matveev, M. A. Salle, Darboux Transfomations and Solitons, (Springer, Berlin, 1991). 

P. Bona, "Quantum mechanics with mean-field backgrounds", Comenius University preprint No. Ph.10-91 (1991). 
T. F. Jordan, Ann. Phys. (NY) 225, 83 (1993). 

I. Bialynicki-Birula and P. J. Morrison, Phys. Lett. A 158, 453 (1991); P. J. Morrison, Rev. Mod. Phys. 70, 467 (1998). 
M. Czachor, Phys. Lett. A 225, 1 (1997). 

M. Czachor and M. Marciniak, Phys. L ett. A 239, 353 (1998). 
M. Czachor, e-print ^uant-ph/9711054; Int. J. Theor. Phys. — to be published. 
M. Czachor, Phys. Rev. A 57, 4122 (1998). 
M. Czachor and M. Kuna, Phys. Rev. A 58, (July, 1998). 
Z. Daroczy, Inf. Control 16, 36 (1970). 
C. Tsallis, J. Stat. Phys. 52, 479 (1988). 
Y. Nambu, Phys. Rev. D, 7, 2405 (1973). 

M. Garcia Sucre and A. Kalnay, Int. J. Theor. Phys. 12, 149 (1975). 
L. Takhtajan, Comm. Math. Phys. 160, 295 (1994). 

J. A. de Azcarraga, A. M. Perelomov, and J. C. Perez Bueno, J. Phys. A: Math. Gen. 29, L151 (1996); 29, 7993 (1996). 
J. A. de Azcarraga and J. C. Perez Bueno, Comm. Math. Phys. 184, 669 (1997). 
J. C. Perez Bueno, J. Phys. A: Math. Gen. 30, 6509 (1997). 

J. A. de Azcarraga, J M. Izquierdo, and J. C. Perez Bueno, J. Phys. A: Math. Gen. 30, L607 (1997) 
G. Ditto, and M. Flato, Lett. Math. Phys. 39, 107 (1997). 

G. Ditto, M. Flato, D. Sternheimer, and L. Takhtajan, Comm. Math. Phys. 183, 1 (1997). 
P. Gautheron, Lett. Math. Phys. 39, 107 (1997). 

R. Ibafiez, M. de Leon, J. C. Marrero, and D. Martin de Diego, J. Math. Phys. 38, 2332 (1997). 
I. Kanatchikov, Rep. Math. Phys. 40, 225 (1997). 
M. Czachor, Phys. Rev. A 53, 1310 (1996). 

P. Grochowski, W. Kaniowski, and B. Mielnik, unpublished (1998). 

The master equation can be easily derived if one keeps in mind that, by definition of the p-inverse, (pxp)~ 1 = p(pxp)~ 1 p 
and 

d(pxp)~ 1 = — (pxp)~ 1 d(pxp)(pxp)~ . 



13 



